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Abstract 

We prove an analog of the Caffarelli-Kohn-Nirenberg theorem for weak solutions of 
a system of PDE that model a viscoelastic fluid in the presence of an energy damping 
mechanism. The system was recently introduced as a possible method of establishing 
the global in time existence of weak solutions of the well known Oldroyd system. 
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1 Introduction 



The Oldroyd model for an incompressible, viscoelastic fluid is governed by the following 
system of equations 

dtu + {u ■ V)m = Au-Vp + V ■ FF* 

dtF + {u-V)F = VuF . (1.1) 

V ■u = 

Informally, we refer to f 1 1.1 1) as the viscoelastic Navier-Stokes equations. This system of PDE 
is always assumed to be satisfied in an open subset of spacetime M.^ x M. At a point (x, t) 
in spacetime, u = u{x,t) G M"^ represents the fluid's velocity, p = p{x,t) G M represents the 
fluid's pressure, and F = F{x,t) G R^^^ represents the local deformation of the fluid. The 
associated energy law for any smooth solution {u,p, F) on x (0, oo) that vanishes rapidly 
enough as |x| — )■ oo is 

d f f |n(a;,t)P \F(x,t)\'^] , /" . m2 , 

for t > 0. 

Solutions of initial value problems associated to (11. ip have been studied extensively. For 
instance, the short time existence of a smooth solution and the global existence of a smooth 
solution that is initially small (in an appropriate norm) has been established in various 
settings |ll|5]. However, it is not known if solutions with smooth initial and boundary data 
develop singularities or even if some meaningful type of weak solutions exist globally in time. 

In pursuing the former problem, the authors of [5] introduced the following system as a 
way of approximating solutions of (11. ip 

dtU + {u ■ V)m = Au-Vp + V ■ FF^ 

dtF +{u-V)F = fiAF + VuF (1.2) 

V ■ n = 

for a parameter /i > 0. The associated energy law for any smooth solution {u,p,F) of (11.20 
on R^ X (0, 00), that vanishes rapidly enough as |x| — ?■ 00, is 

for t > 0. Therefore, the presence of /i acts to create energy dissipation, and so we interpret ^ 
as a damping parameter and call (II. 2p the viscoelastic Navier-Stokes equations with damping. 

It is not difficult to establish the existence of a global in time weak solution of (ll.2p 
analogous to that of the Leray-Hopf solutions of the incompressible Navier-Stokes equations. 
Unfortunately, standard weak convergence methods do not allow one to pass to the limit as 



2 



/i — )■ 0"'' to generate weak solutions of (11. ip : see the end of section 2 in reference [5] for more 
this. Nevertheless, the system f ll.2p is itself of interest and is the topic of study in our work. 

A fundamental simplification that we will make in our analysis is that we consider solu- 
tions {u,p, F) of (11. 2p that additionally satisfy the equation 

V-F* = 0. (1.3) 

That is, a standing assumption that we shall make is that the deformation F has diver- 
gence free columns. This assumption is motivated from taking the divergence of the second 
equation in (II. 2p which yields the following transport equation 

dt{V ■ F*) + (n • V)(V ■ F*) = /iA(V ■ F*). 

The above PDE formally implies that if (II. 3p holds at some instance of time, then it will hold 
at all later times. Therefore, we believe that our results below will be pertain to solutions 
of that initially satisfy ffO]) . 

As our results do not change qualitatively as > is varied, we set 

/i = 1 

in our analysis of weak solutions of (II. 2p . Our main result is the analog for (II. 2p of the 
Caffarelli-Kohn-Nirenberg theorem for the incompressible Navier-Stokes equations [2]. The 
statement of this theorem involves the concept of a weak solution, which we will define in 
the next section, and the concept of the singular set of a solution {u,p, F). The singular set 
corresponding to a solution {u,p,F) is defined as the set of points {x,t) in the domain of 
{u,p, F) for which either u of F is not Holder continuous in any neighborhood of {x, t). Any 
point not belonging to the singular set of {u,p, F) is a regular point. 

Theorem 1.1. There is a universal constant e > such that if {u,p,F) is a weak solution 
of (O and zf 

limsupi// {\Vu\^ +\VF\^}dyds <e, (1.4) 

r->0+ 1^ J JQr{x,t) 

then u and F are Holder continuous on some neighborhood of{x,t). In particular, {x,t) ^ S. 

In the limit above (II. 4p . and in this work, 

Qr{x, t) := Br{x) X {t- r^/2, t + r'^/2) C x R 

is a parabolic cylinder of radius r centered at {x,t). Using standard covering arguments (see 
in particular section 6 of [2]), we have the following corollary. 

Corollary 1.2. Assume that {u,p,F) be a suitable weak solution of (ll.2p on an open subset 
o/ R'^ X R and let S be the singular set of u and F. Then V^{S) = 0, where denotes 
one- dimensional parabolic Hausdorff measure on M.^ x R. 
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In proving Proposition ll.H we employ the blow- up/compactness method introduced by 
Lin [3] that lead to a relatively simple proof of the Caffarelli-Kohn-Nirenberg theorem. We 
also borrow many of the ideas from the very clear account on these topics given by Ladyzhen- 
skaya and Seregin |6]. The organization of this paper is as follows. In section |2l we define 
suitable weak solutions and establish an important compactness result for these solutions; 
a fundamental corollary of this compactness is a "decay" or "blow-up" lemma. In section 
El we use the decay lemma to deduce a local condition that if satisfied implies solutions are 
locally Holder continuous. Finally, in section HJ we show that (11. 4p implies that our local 
regularity condition is satisfied, which in turn furnishes a proof Theorem 11.11 

2 Weak solutions 

As mentioned above, we set = 1 in equation (II. 2p and study the system of PDE 

dfU + {u ■ V)m = Au-Vp + V ■ FF* 

dtF + {u-V)F = AF + VuF . (2.1) 

V ■ n = 

Sometimes it will be beneficial for us to write equation (12. ip in terms of the columns of the 
matrix valued mapping F. Setting Fj := Fej for j = 1,2,3, we have from the assumption 

3 

V-FF* = J](F,-V)Ffc. 

fc=i 

In particular, (12. ip and (11.30 can be rewritten together as 

dtu + (n • V)n = Au-Vp + ELi(^fc " V)Ffe 
dtF, + {u-V)Fj = AF, + {Fj-V)u, j = 1,2,3 . (2.2) 
V ■ n = V • = 

As Theorem 11.11 is local, we only consider solutions on the unit cylinder Qi := Qi{0,0). 
Before pursuing the analysis of solutions, let us make some basic observations that will mo- 
tivate the definition of suitable weak solutions and other ideas to follow. 

Scale invariance. If {u,p,F) is a solution of (12. ip on Qr{xQ,tQ), then {u^,p'^,F^) is a 
solution on Qr/x{0, 0) for A > 0, where 

u^{x, t) = Xu{xo + Ax, to + A^t) 

u^{x, t) = X^p{xq + Ax, to + A^t) . (2.3) 
F^(x, t) = XF{xo + Ax, to + A^t) 
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Local energy identity. If {u,p,F) is a smooth solution of (12. ip on Qi and G C^{Qi] 
then 



Bi 



n|2 




2 


2 






2 ^ 


2 



for t G (-1/2,1/2). 



A spacetime bound on n, F. From the local energy identity above, we expect 

solutions {u,p,F) to satisfy 



sup / \u\^ + \F\'^dx+ 1 1 \Vu\^ +\VF\^dxdt<C. 

-l/2<i<l/2 Jbi 



An application of the interpolation estimate 



for V G H'^{^l;R^), where 
provides the bound 



3 1 

a = , 2<r<6, 

r 2' - - ' 



1 10/3 



\ur'- + \F\^^/^dxdt<Ci. 



Qi 



Here Ci only depends on the constant C in 
A spacetime L^/^ bound on p. Taking the divergence of equation (12. ip gives 



- Aj9 = V 



V)u-J2iFk-V)Fk 



k=l 



Note that 

|(m ■ V)M|ii5/i4(B^) < |VM|L2(Bi)|M|i30/i3(Bi)- 

By the interpolation inequality (I2.5p with r = 30/13 and a = 4/5 



|«Il3o/i3(bo < c||Vn|^/^(^^^ + l| . 
With dMl) and ([23]), we have 

|(«-V)«|^S/14(^^)<C{|V«|i.(^^) + l} 
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(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2.^ 



(2.9) 



and likewise 

■ V)F,|^5/,,(^^^ < C {|VF,|i.(^^) + 1} , 3 = 1,2,3. 

Since, 

'15\* 5 



14/ 3' 

we also have by the Sobolev inequality and basic PDE estimates following from equation 



[ [ \p\'/'dxdt <C f'^ \Wp\'l' dt 



<Cjj |^|Vn|2 + J]|VFjf + lj dxcit (2.10) 

for any 9 G (0, 1). Thus p E ivf/f (Qi). 

The above observations motivate the following definition of weak solutions of the system 
(12. ip . This definition is of course also consistent with the notion of (suitable) weak solutions 
presented in the original work by Caffarelli, Kohn and Nirenberg [2]. 

Definition 2.1. {u,p,F) is a weak solution of ( 12. ip on Qi provided 
(i) u,Fe L\Qi), and p E I^"\Qi), 

{ii) equation (12. ip holds in the sense of distributions on Qi, and 
{in) for each E C,°°(Qi), with > 

/ ^f^ + ^V^+r [ <P{\S/u\' + \\/F\')dxds< f /■ /(0, + A0)f^ 
JBix{t} V ^ ^ / J-1/2JB1 J-1/2JB1 L V ^ 

12 I 1712 
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-u(g)V0-FF*+ + n- V0|rfxds (2.11) 

te [-1/2,1/2]. 

Our first order of business is to establish that a bounded sequence of weak solutions has 
a convergent subsequence whose limit is again a weak solution. We will not use this result 
directly altough the ideas that go into proving this result will be essential to our proof of 
Theorem 11.11 To this end, we start by quoting a relatively standard compactness lemma. 
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Lemma 2.2. Let Xq, Xi, X2 denote Banach spaces, with Xq and X2 reflexive, that satisfy 

Xq C Xi C X2. 

Also suppose that the embedding of Xq into Xi is compact and the embedding of Xi into 
X2 is continuous. Let p,q G (l,oo) and assume that {uk)k£N G -^'^([0, 1]; ^0) ^-5 a bounded 
sequence such that each Uk has a weak derivative u'^ and the sequence 

{u'k)kmeL\%l]-X2) 

is also bounded. Then there is a subsequence of Uk converging strongly in L^([0, l],Xi). 

Remark 2.3. We can replace [0, 1] in the statement above with any compact interval of M. 

We refer the reader to Theorem 2.1, section III of [8] for a detailed proof of the above 
lemma. Its primary application is the following theorem. 

Theorem 2.4. Suppose that {u^,p'^,F^)k<£n is a sequence of weak solutions of (12.11) on Qi 

satisfying 

{\u''f + \p''\^/^ + \F''\^} dxdt < C (2.12) 



for all k > 1 and some constant C > 0. Then there is a subsequence [u^^ ,p^^ , F^^)j(z'M and 
{u, P, F), such that 

u^^ -^u in Ll^^{Qi) 
pk, L^/^iQi) 
F^.^F in LUQ,) 
for 1 <q < 10/3. Moreover, {u,P,F) is a weak solution of (12. ip on Qi. 

Remark 2.5. It is straightforward to adapt the proof below to build globally defined weak 
solutions of a large class of initial value problems associated with (12. ip . One may use a 
Galerkin-type approximation, for example. 

Proof. First, we select the weak limit P of the sequence {p^)k&n and without any loss of 
generality assume p^ P in L^^'^{Qi) as /c — )■ cx3. Next, we observe that by (12. 2p . we have 
that 

dtu'' = V ■ {-u'' ®u'' + Vu^ - p^L + F^{F^Y) 

and 

dtF^ = V ■ ® m'" + ® Ff + VF/) , J = 1 , 2, 3 

belong to the space 

L^I\{-l/2M2)-W-^'\B,)). 

By (12.121) . dtu'^ and dtF^ are in fact bounded in this space for all > 1. 

For a fixed 9 G (0, 1), we note that as {u^^p^^ F^) satisfies the local energy estimate (12. lip 

iVn^f + \VF^\^dxdt < C (2.13) 



for some C = C{6). These observations lead to the choice of exponents 

p := 2 and q := 3/2 

and Banach spaces 

^ H\Bo) 

By our observations above, the hypotheses of the previous lemma are satisfied with 
the sequences {u'')k£N and {F'')k(^n (and exponents and Banach spaces as indicated above). 
Hence, some subsequence of (u^)^^^, (-F'^)^^^ converge in L'^{Qg) to some ug,Fg. By the 
bound (12.131) . u^.F'^ are bounded in L^^/'^^Qg) (recall the estimate (12. 6p ). Consequently, 
the interpolation of Lebesgue spaces implies {u^^F^) — >■ {ug,Fg) as /c — )■ in L'^{Qg) for 
1 < g < 10/3. As ^ G (0, 1) was arbitrary, we can employ a routine diagonalization argument 
to construct a u, F G L\^[^{Qi) such that a subsequence of and F^ converge respectively 
to u and F in Ll^^{Q) for 1 < g < 10/3. It is now immediate to pass to the limit as A; — )■ oo 
and show that [u, P, F) is a solution of (12. ip in the distributional sense and that (12. lip holds. 
Therefore, {u, P, F) is a weak solution of (12. ip on Qi. □ 

Another application of the compactness lemma is the following "blow-up" or "decay" 
lemma; these names come from the fact that the lemma's proof uses a rescaling and blow-up 
argument, while the conclusion involves a type of decay. This is arguably the most important 
step in the proof of Theorem 11.11 These ideas originated with the groundbreaking work of 
Lin [3]. However, the specific approach we use here follows closely the reinterpretation by 
Seregin and Ladyzhenskaya [6]. 

Lemma 2.6. Let {u,p,F) be a weak solution on Qi and set 

E{x,t,r) := (^j-J]u - Ugfdyds^ +r (^j-j\p - PBAt)f^^dyds^ ^(/^^ ~ FQrfdyds^ 
for Qr = Qr{x,t) and Br = Br{x), where 

UQr ■= 4-ru{.y^s)dyds and PB^is) := +p{y, s)dy. 

For each 

0<e <l/2 and M > 
there are positive numbers ei,i?i,Ci such that if 

{i) Qr{x,t) C Qi, r < Ri 

< {ii) \uQ^{x,t)\ + \FQ,{x,t)\ < M , 
Hi) E(x, t, r) < ei 
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then 

E{x,t,er) < Cie^^^E{x,t,r). 
Moreover, c\ can he chosen independent of 6. 

Remark 2.7. In particular we can arrange c\(f'l^ to be small, so this is a type of decay lemma. 

Proof. 1. We argue by contradiction. Suppose the statement of the lemma is false. Then 
there are sequences 

Cfe ^ 
rfc ^ 

^ Cfc = Ci (which will be chosen sufficiently large below) 



with 



and 
We set 



{%) Qrk{Xk,tk) C Qi, 
' (^0 \uQr,{x„t,)\ + \FQ,.^ix„t,)\ < Mo 
^{iii) E{xk,tk,rk) = 
E{xk,tk,Oork) > ci6l^'^E{xk,tk,rk). 



(2.14) 



(2.15) 



\''{y,s) ■=e^^ {u{xk + rky,tk + rls) - UQ^.^(^^^^t,)) 

< p^{y,s):=rkel^ {p{xk + rky,tk + rls)-pBr^{xk){tk + rls)) , (y, s) G Qi 



f'iy^s) := e;.^ (F(xfc + rfc|/,tfc + r2s) -Fq^j,,,^,)) 

and observe that {u'',p'^,F^) is a weak solution (defined analogously as in Definition 12.11) of 
the PDE 



dtu'' + ((a^ + ekU^) ■ V)m^ = Au^ - Vp^ + EjLi((^^ + Cfc^/) " V)F/ 
dtF^^ + ((a^ + eku'') ■ V)F/ = AFf + ((6^^ + e^F/) ■ V)m'= 



V ■ = V ■ F^' = 



{y,s) e Qi. 

(2.16) 



for j = 1, 2, 3. Here 



a'" ■■= UQ,Jx„t,) and 6j := (F,)q,Jx„ 
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are bounded sequences in for each j = 1,2,3 by fl2.14p . Moreover, the sequence {u^, p'', F^) 
satisfy the integral bounds 



u 



fe|3 



1/3 



Qi 



2/3 



fc|3 



1/3 



efc 



and the generahzed energy inequahty 

k\2 |77fc|2 



Bix{s} 



dx+ r [ 0(ivu'=|2+ivFY)d2;(is< r / 



-1/2 7Bi 



1/2 7Bi 
^fc|2 |^fc|2 



(2.17) 



Ufc|2 |^fc|2 
+ 



+ 



2 2 



(2.18) 



dxds 



for G C,°°(S), (/) > and s G (-1/2, 1/2). 

2. Arguing as we did in the proof of the Theorem 12.41 one checks that the sequence 
{u^)k&i and (F*^)fcgN satisfy the hypotheses of the Lemma [2l2l Moreover, we conclude that 
there is m G Ll^^{Qi), p G Lfl^{Qi), and F G Lf^^{Qi) and subsequences of {u'')kmi {p'')keN, 
and (-F*^)fceN (again labeled u^,p^,F^) such that 

^ M in Lf^^iQi) 
p^ ^ p'm. L^/^(Qi) 
F'^^Fin FL(gi) 

as A; — oo. From f l2.14p . we may assume without any loss of generality that that 

-.k 



a ^ a 
6? ^ bi 



J = 1,2,3 



in M^, as — )■ cxd. 

It is follows from this convergence that {u,p, F) are weak solutions of the following linear 
PDE 

dtu + (a ■ V)n = Au-Vp + Ei=i(&i ■ V)Fj 
^tF,- + (a ■ V)Fj = AFj + (6^- ■ V)u (2.19) 



V ■ M = V ■ F, = 



and satisfy the energy inequality 
n|2 |F|2 



Bix{s} 



2 ~^ 2 



^■^)dx+ r [ (i){\vu\^+\\7F\^)dxds< r [ 

^ J J-I/2JB1 J -1/2 J B 



+ < 



1/2 JBi 

2 I t:ii2 



^ + ^H<^* + A<^) 



J](n • Fj)bj + [LL + i^) a + pu \ ■V<P 

i=i 



dxds 



(2.20) 
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for e C~(S), > and s e (-1/2, 1/2). 

3. We claim that there is a constant Ci = Ci{Mq) such that 

\u{x,t)-u{y,s)\ < Ci{\x - y\ + \t - s\^/^) 
\F{x,t)-F{y,s)\<C,{\x-y\ + \t-s\), 

for {x,t), {y, s) G Qi/2- In particular, we are asserting that u is Holder continuous and F is 
Lipschitz continuous on Qi/2- 

To see this, we first take the divergence of the first equation in (12.191) to get for almost 
every t G (-1/2, 1/2) 

-Ap(x,t) = 0, X e Bi. 

Thus p{t) G C°°{Bi) for almost every t G (—1/2,1/2). Recall the following estimate for 
harmonic functions: —Ah = 0, in Bi C M", then for each multiindex a, q > 1, and 6 G (0, 1), 
there is a C{q, a) such that 

|5"/i|L-(Be) < -^Y^^fy^l^linsi) (2-22) 

(see Theorem 7, page 29 of [1]). As p G L^^'^{Qi), we conclude that 

9:pGL=^/^((-l/2,l/2),L-(i?0) 

for each multiindex a. 

Next, we take the curl of the first two equations in (I2.19P to arrive at 



where of course 



dtw + (a ■ V)w = Aw + Ej=i(^j ■ V)rj 
9jrj + (a ■ V)rj = Arj + {bj ■ \/)w 



w := W X u and rj := V x Fj. 



(2.23) 



From the local energy estimate for {u,p,F) (12.201) . we have w,rj G L^^^{Qi). As 

dtW-Aw = V- (Ei=i rj (^bj-w(E)a^ e L\{-l/2, 1/2), H'\B^)) 
dtVj - Arj = V -{w^bj- rj ® a) G L^ii-l/2, 1/2), H~\Bi)) 

It follows from standard energy estimates for the parabolic system (I2.23P that 

w, r, G L°°((-l/2, 1/2), L\B,)) H L\{-1/2, 1/2), H\B,)). 

As (12.231) is linear, and in particular each of the derivatives d^w, d^rj also satisfy (I2.23p . we 
conclude by induction that for each multiindex a 

a>, d^r, G L-((-l/2, 1/2), L\B^)) n ^^((-1/2, 1/2), H\B,)). 



11 



Therefore, each space variable of both w and vj is locally bounded on Qi. 
Fom the condition V ■ m = 0, we have 

- Au = V xw. (2.24) 

In particular, the Biot-Savart law reads for for each compact G G Bi, 

u{x,t)= / V(f){x -y) X w{y,t)dy + A{x,t), (x, t) G G x (-1/2, 1/2) 



JG 

where G 3 x ^ A{x,t) is harmonic for almost each t G (—1/2, 1/2) (see Lemma 2 of [7] 
for a proof). From the local energy estimate for the solution {u,p,F) (12.201) . we know that 
the LP'{G) norm of u is bounded independently of t belonging to compact subintervals of 
(—1/2, 1/2). From our estimates on we conclude that the LF'{G) norm oi A is also bounded 
independently of such t. As A is harmonic, it is immediate from the mean value property 
that A is locally bounded on F x I, for each interval I C (—1/2, 1/2). Hence u is locally 
bounded on Qi; differentiating (I2.24p and recalling our estimates on w, we also see that d^u 
is locally bounded on Qi. 

Arguing in the same manner, we have that d^Fj is locally bounded on Qi for each 
multiindex a and j = 1, 2, 3. It is now immediate from equation (12.191) that 

3 

d,u = -(a ■ V)n + An -Vp + J^i^j ' V)F, G L?/,'((-l/2, 1/2), 1^(8,)) 

i=i 

and 

dtF, = -{a- V)F, + AF, + (6, ■ V)n G ^^^((-1/2, 1/2), LZiB^)) 

from which the claimed estimates (12.211) readily follow. 

4. Therefore we have that there is constant G2 = G2{Mq) such that 

u - UQ,^^ \Hyd^ + (^jj-W - Fqoo I'dyd^^ < C^^T ■ 
By the convergence in L\^Jy(^^ of [u^ , F^) to (u, F) as A; — )■ 00, 

^ - n^^^ I'dyds^ + (^fj^\F' - F^^^ fdyds^ < IG^bT (2-25) 

for all /c G N sufficiently large. 

Let us now derive an estimate on the average of the p^ . Direct computation from (12.161) 
gives 

3 

-A/ = efcV ■ \{u^ ■ V)n^ - J^iP' " V)i^1, x e B, 
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for for almost every s G (—1/2,1/2). For each s G (—1/2,1/2) that this equation holds 
weakly on Bi, let g^'i^s) denote the unique solution of the PDE 

-Ag'^ = e,V ■ [{u'^is) ■ V)u^{s) - ELi(^/(s) " V)F/(s)], x G ^3/4 



g'' = 0, X G dB: 



3/4 



and set 

,= pk^s)-g'{s). 

Clearly, h'^^s) is harmonic on -B3/4. Virtually the same argument that lead us to the bound 
fl2.10p provides the estimate 

/ / Ig'f'dyds < Cek 

J JQi/2 

where C depends (say) only on an upper bound for the integrals / Jg^^ (| Vu'^p + | V-F'^p) dyds. 

The local energy estimate for the solution {u^,p^,F'^) fl2.18p combined with the Li^{Qi) 
bounds f l2.17p for (m'^, F^) assures us that these integrals are uniformly bounded above. 
Therefore, by these remarks and the estimate f l2.22p 

''W^^-^n ^-'iWf-<n HW!f-''-n 



^ ^ [IL [kl°Us^''''''-^'f 



2/3 



X-'-^Ws/s 
2 

7/3 ' """0 
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Consequently, for k sufficiently large 

\ 2/3 

Ip'-PBsf^'dyds] <C,el/' (2.26) 

So / 

for some universal constant C3 (since for all large k we have < 6q). Combining f l2.25p and 
f l2.26p . it is readily checked that a contradiction to fl2.15p is obtained, for all k sufficiently 
large, by choosing 

Ci(Mo) := 2(2C2(Mo) + C3), Mq > 0. 
Moreover, our choice Ci is independent of 6. □ 



3 A local criterion for Holder continuity 

With Lemma [2.61 in hand, we are now in position iterate its conclusion which is the major 
step in establishing a local criterion for Holder continuity for weak solutions of the system of 
PDE f l2.ip . As a corollary of this Holder continuity criterion, we prove Theorem 13.51 which 
is reminiscent of Serrin's higher regularity result for weak solutions of the incompressible 
Navier-Stokes equation [7]. 

Lemma 3.1. Assume {u,p,F) is weak solution on Qi and fix 

0<9 < 1/2, M > 0, < /3 < 2/3 (3.1) 

so that 

Ci(M)el/3-/3/2 < ^ 

There is €2 > such that if 

\i) Qr{x,t) cQu r < Ri 

' (^0 \UQr{x,t) \ + \FQ,{x,t) \ < \M 

(Hi) E{x, t, r) < €2 

then for any k = 0,1,2,3, .. . 

(0 \FQg,^(x,t)\ + \uQ^,^i^,t)\ < M 

< (ii) E{x,t,e^r) < ei 

^ {Hi) E{x, t, e^+^r) < ^(^=+1)^(1 - e^l^-PI^)-^E{x, t, r) 
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1/3 / /. /. \ 1/3 

|3^.,^„\ , (ll\p^ ^ 13, 



I(x,t,r) 



Proof. 1. First set 

I{x,t,r):= i-U\u-UQfdyds] + [-f-f-lF - Fg^l'-'dyds 

\J JQr 

for Qr = Qr{x,t), and notice that 
The above estimate follows directly from Holder's inequality, and using the triangle inequality 

^ k-l 

We define 



j=0 



r{x,t) I 



e, := (1 - ^1/3-/3/2) min jiei, ^^/^(i _ 

as our candidate for e2 described in the statement of the lemma. We argue by induction 
below. 

2. We establish the claim for k = 0. By assumption and our choice of 62 
E{x,t,r)<e2<{l-e'/'~^/'fj<'j<e,. 
It is also easy to verify that the conditions of the Lemma 12.61 are satisfied, and thus 



E{x,t,er) < Cie^^^E{x,t,r) 

< c,e'/^-^/'e'/'+''/'E{x,t,r) 

< l-9'/'+''/^E{x,t,r) 

< 9^(l-e'/''-^/')-'E{x,t,r) 



3. Assume for s = 0, 1, . . . , /c 



as ^i/3+/3/2 < _ 0i/3~p/2yi ^^j^i^j^ happens if and only if ^V3"/3/2 < 3). 

(,'i)s \uQgs,{x,t) \ + \FQgs,ix,t)\ < M 
< E{x,t,e'r) < ei 

^ {ni)s E{x, t, P+V) < - ^1/3-/3/2) -1^^^,^ ^) 

We show that these assumptions imply that the above bounds hold in turn for s = k + 1. 
For simplicity, we suppress the (a;, t) dependence in our arguments below. 
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For {i)k+i- 

k 



< 



3=0 



< J-j(l-»'/3-«.)-.E(r)X:C» + f 



< 



Q5/3 

j=0 

1 E{r) M 

^5/3 (1_ ^1/3-/3/2)^1 + Y 

< ^ + ^ 

- ^5/3(1 _ ^1/3-/3/2)^1 ^ 2 

M M 

^ T + T 

= M. 



For By assumption {ni)k, we have 



< ^(fc+l)/3(i_ ^1/3-/3/2)-!^^ 

2 ' 

< ei- 

For {in)k+i. Observe that the hypotheses of the Lemma (12. 6p are satisfied with r 9^r. 
Moreover, using (13. ip 

< Ci^2/3^(^-+l)'^(l - ^l/3-/3/2)-l^(^) 

^ ^^^l/3-/3/2^1/3+/3/2^(fc+l)/3(^l _ ^l/3-/3/2)-l ^(^^) 

< 1 . 5)l/3+/3/2^(fc+l)/3(^i _ ^l/3-/3/2)-l ^(^^) 

< ^(fc+2)/3(i_^l/3-/3/2)-l^(^)^ 

The last inequahty follows since /3/2 + 1/3 > /3 and 6 G (0,1) which trivially implies 

5)l/3+/3/2 < ^/3^ 

□ 

The main use of the decay and iteration lemmas is the following proposition. We omit 
the proof as it follows from the above lemma and relatively standard manipulations. See 
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Proposition 2.8 in [6] for a related result, concerning the incompressible Navier-Stokes, and 
its proof, which is easily adapted to our framework. 

Proposition 3.2. Assume that {u,p,F) is a weak solution of (12. ip . There are universal 
constants e^, R2 such that if 

Qr = Qr{x,t) C Qi, r < R2 
E{x,t,r) < es 

then u and F are Holder continuous in some neighborhood of {x,t). Here 

E{x,t,r) := fZ/l^I'V^' + r f//bp/^y^% f//|F|=^y^'. 

The scaling invariance properties of the system (12. ip imply the following improvement of 
Proposition fl3.2p F1 This simple observation is as important as any we make in this work. 

Corollary 3.3. Assume {u,p, F) is a weak solution on Qi and let €3, R2 be as in Proposition 
\3.2[ Further suppose that 

(Qr{x,t)cQu r < i?2 
\^E{x,t,r) <e3 

Then {x,t) is a regular point for {u,p,F). 
Proof. Set A = r/i?2 and define 

u^{y, s) := Xu{x + Xy,t + A^s) 

u\y,s) := X'^p{x + Xy,t + X^s) ,{y,s) e Qr,{0,0). 
F^{y, s) := XF{x + Xy,t + A^s) 

It is immediate that is a weak solution on Qji2{0,0), and direct computation 

yields 

E(0,0,i?2;M\p\F^) = XE{t,x,r;u,p,F) 

r 

= —E{t,x,r;u,p,F) 

< 63- 

By Proposition 13.21 and F^ are Holder continuous in a neighborhood of (0,0). Conse- 
quently, u and F are Holder continuous in a neighborhood of {x,t). □ 

1 Recall dlH). 
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Corollary 3.4. Assume {u,p, F) is a weak solution on Qi and let 63, R2 be as in Proposition 

\M If 

liminf rE{x,t,r) < -e^R2 

r— >-0+ 2 

then {x,t) is a regular point point for {u,p,F). In particular, there is €4 > such that the 
same conclusion holds provided 

liminf 4 /"/" \\u\^ + \p\^^^ + \Ff] dyds < 64. (3.2) 

Proof. If 

liminf r-E'(x, t, r) < -e^R2, 

r-->0+ 2 

then info<r-<<5 ^, < 1^3-^2 for any 5 > 0. Moreover, there is r < i?2 such that 

■^E{x,t,r) < es and Qr{x,t) C Qi. The previous corollary then applies. The second 

assertion follows immediately from the first. □ 

The following claim is in the spirit of Serrin's regularity criterion for weak solutions of the 
incompressible Navier-Stokes equations. Note, however, that our result requires integrability 
of the pressure as well as integrability of the velocity and the deformation. 

Theorem 3.5. Assume that {u,p,F) is a weak solution on Qi such that 

"1/2 / p \ s'/s 



1/ 

for some s' > s satisfying 



/•I/2 / r \«/« 

/ / + \p\'^^ + \F\'} dx ) dt<oo. 

i-1/2 \Jbi J 



(3.3) 



3 2 ^ 

- + -<1 

s s' 



Then u, F is regular on Qi. 

Proof. For each point {x,t) G Qi, the bound ( 13. 3p implies 

lim^ [[ \\u\^ + \pf^ + \F\ndyds = 0, 
^^or^ JjQ,{x,t)' 

which in turn implies (13. 2p . Therefore, the claim follows from the previous corollary. □ 

4 "A-B-C-D" estimates 

In the previous section, we deduced that there is an es > such that if 

liminf r£'(a;, t, r) < €5, (4.1) 

r—>-0+ 
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then u and F are both Holder continuous near (x, t)o We claim that the central hypothesis 
of Theorem 11.11 implies the above limit. More precisely, we assert the following fundamental 
proposition. 

Proposition 4.1. Let {u,p,F) be a weak solution on Qi. There is a universal constant 
e > such that inequality (14. ip holds whenever 

limsup - // {|Vnp + \VF\^] dyds < e. 

r-i-0+ ^ J JQr{x,t) 

In view of Corollary 13. 4[ a proof of Proposition 14. II establishes Theorem ll.il To this end, 
we shall need three estimates involving the following integral quantities. For (x, t) G Qi, and 
r > so small that Qr{x,t) C Qi, we define 

A{x,t,r) := sup|i_,j<^2/2i/5^(^.){|u(y,s)p + |F(y,s)|2}c??/ 
B{x, t, r) := i //^^(^^^^ {\Vu{y, s)\' + \VF{y, s^} dyds 

< 

C{x,t,r) := ^ JjQ^^^^^^{\u{y,s)\' + \Fiy,s)ndyds 

. D{x,t,r) := ^//^^^^.^^ \p{y,s)\^/^dyds 

Our arguments below are independent of (x, t), so without any loss of generality we establish 
our results for = (0,0). For ease of notation, we also write A{r) := A(0, 0, r), _B(r) : = 

B{Q,0,r),C{r) := C(0,0,r), and D{r) := D{0,0,r) for < r < 1. 

We will need three estimates before undertaking the proof of Proposition 14.11 They are 
very similar to the sequence of Lemmas needed in [3] and as in the previous section we follow 
the path of [6] closely. The statement of the lemmas in fact are nearly identical to the ones 
use to prove the version of Proposition 14.11 in [6], but unfortunately the proofs had to be 
modified. Nevertheless, we would like to emphasize that the work of [3] and [B] served as a 
guiding light for this section. 

Lemma 4.2. There is a universal constant c > such that 
for < r < p. 

Proof. This assertion follows directly from following the well known inequality established 
in [3] (Lemma 2.1). There is a universal constant c > such that 

Coir) < c| (P) + (^)'^f (P)50^'(P)| > < r < p < 1 

^£5 := ^£3-^2 
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for all 
where 



V e H\Q,) n L~((-l/2, 1/2); L\B^)), 



Ao(r) 
Coir) 



sup|fi<^2/2^/5^ \v{y,s)\'^dy 
T^IIo^ \v{y,s)\'^dyds 



□ 



Lemma 4.3. There is a universal constant c > such that 

A{r/2) + B{r/2) < c {^^/^(r) + C^/''{r)D^/%r) + A^^\r)B^^\r)C^/%r)} 
for < r < 1 . 

Remark 4.4. Recall the energy inequality fl2.1ip : G C^{Qi),(j) > 
m|2 



Bix{t} 



-1/2 JBi 

-F'u ■ F'V(f) + 



0(|Vm|2 + |VFp) dxds < 



1/2 JBi 



A0) 



t e [-1/2,1/2]. As 



V ■ M = and V ■ F* = 0, 
we can replace the right hand side of the above inequality by 

for any a G 1/2, 1/2; M.^) and 6 G Ll^^{~l/2, 1/2). In a crucial step below, we will use 



/' /{ 




1-1/2 J Bi I 





'a(t) := [F%]BXt) = fBr^\x,t)u{x,t)dx 
b{t) := iHlit) + l\F\lit) = 



Hx,t)\^ , \F{x,t)\^ 
2 2 



(4.2) 



Proof. Choose </> G C^{Qi), nonnegative, such that = 1 on Q1/2, and set 

' X t 



(j)r{x,t) := (p 



' 9 I ' 



(X,t) G Qr 



Clearly 4>r G C;?°(g,.) and 

'|V</)r| < c/r 
|ai</)^ + A(j)r\ < c/r'^ 

for a constant c > independent of r G (0, 1). 
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Using 0r as the test function in the energy inequahty (12. lip and choosing the mappings 
a and b as in (14.21) gives 

A(r/2) + i?(r/2)<c(l // iH' + \F\') + ^ [ [ 11^^-1^11 l|n| 



1 



|3 



Qr 

r 

Qr 

1/3 / ^ /. N 2/3" 



^ \\F\' - \F\IM + II \\F'u\ - [F'uU.\\F\ 



\pf') ■ (4.3) 



We now proceed to estimate each integral in the inequahty above (14. 3 p in terms of A, B,C,D. 
First note that 



Uf + im < ( /■/■ M'+iFi'V"iQ.i'" 



2/3 

<Cr5/3-3( II (|^|3^|^|3)- 



= cC(r)2/3, 

for some universal constant c > 0. Next, observe 
We also have, by employing Holder's and Poincare's inequality. 
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N^-HIJH< / ( / IH^-Hl.r ) ( / \un dt 



-r2/2 VJBr / \JBr / 

< I / |m||Vm| / \uf) dt 

rrV2 / f \ / /• \ 1/2 / /• \ 1/3 

< / / i«n / / i^n 

-r2/2 \JBr J \JBr J XJBr J 

rrV2 / c / r \ 1/3 



<cr'/^A{rf" / ( / \VuA ( / \uA dt 

sc,.vMM'-(//iHf'7r(/jv„r)"*) 



1/2 

2 



-r2/2 V^Br / / 

< Crl/2+2/3^(^)l/2cl/3(^) /"/■/■ I V«|2^ //2 

< Crl+2/3^(^)l/2^1/3(^) I 

< Crl+2/3+l/2^(^)l/2cl/3(^)5l/2(^) 

< cr2+i/6A(r)i/2Ci/3(r)5i/2(r) 

< cr2A(r)i/2Ci/3(r)fii/2(r), 

for some universal constant c independent of r G (0, 1]. Arguing in virtually the same manner 
shows 

[[ \\F\'- \F\l.\\u\ + \\F'u\ - [F'u]Bm < cA{rf''C''\r)B'l\r). 

Substituting all of these bounds into inequality f l4.3p gives the desired estimate 
A{r/2) + B{r-/2) < c {C2/3(r) + C^/^^s/s^^) ^ A{ry/^C^/'^{r)B^/^{r)} . 

□ 

Lemma 4.5. There is a universal constant c > such that 



Dir)<c!^^Dip)+[^yA'/\p)B'/\p) 



< r < p. 

Proof. 1. Taking the divergence of the first equation in (12.1 1) gives 

- Aj9 = V ■ [V ■ {uu' - FF')] = V ■ [uu' - [uu']p - {FF' - , (4.4) 
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as (m ■ V)m = V ■ MM*. Here [mm*]^ denotes the average of uu^ on Bp and likewise for [FF^]p. 
Now let Qi = Qi{t) G Wq''^^^{Bp; M^^^) denote the unique weak solution of the PDE 

-AQi = nn* - - (FF* - x G 5^ 

Qi = 0, x G ' 

for (almost every) |t| < 1/2. It is immediate that 

Pi = V-[V-Qi]eL^/^{Bp) 

satisfies (14. 4p . and by the Calderon-Zygmund inequality 

IpiIls/hb,) < |V ■ [V ■ Qi]L:v2(B,) < C{\uu' - K],L3/.(B,) + \FF' - [FF%|i:v2(B,)). (4.5) 
It is also clear that 

P2 ■.= p-pi 

is harmonic in Bp, for almost every \t\ < 1/2. 

2. By the estimate (14.50 and Poincare's inequality, 

2/3 

^ / rill » r I I 7— T I I » — r 7— T I T , 



I \pi\'^'^dx\ <cl {\u\\Vu\ + \F\\VF\}da 



1/2 / X 1/2 

<c| ^ {\u\' + \F\')dx\ I J (|Vnp + |VF|2)dx 

2^ 



Integrating time from — p^/2 to +p^/2, gives 

// bi|'/'<cp3/'A=^/^(p) r^'f/ (|V«p + |VF|2)rfx) 

^^Qp \JBp J 

( \ 

< Cp3/4A3/4 ^ y ^ (|Vlz|2 + I VFHrfxrft j pl/2 



< cp^A''I^B[pfl\ 

It also follows that 



^^Qp >/^Qp >/^Qp 

<cp^{D{p)^A^I'B{^pfl'). 
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3. p2 is harmonic and thus |p2|^''^ is subharmonic for almost every |t| < 1/2. In particular, 

r ^ / \p\^/'^dx 

^ J Br 



is nondecreasing. Consequently, 



1 /■ /■ , ,-:./o r 1 



oil IJ^^I — 9 



for < r < p. Combining the preceding inequalities yields 
D{r) = \ [[ \p\^^^dxdt 



< 

r 



< 
r 



-J I \p,f'dxdt + -^ II \p,fHxdt 

JjQr ^ JJQt 

-II \p,fHxdt + ---^ II \p,\'/'dxdt 

JjQp P P JjQp 

< [P-^ A'/\p)B'/\p) + J- [D{p) + A'/'B{pf/'] 
<^[-D{p)+{^)'A'''B{pf'']. 



as desired. □ 

We are finally ready to prove of Proposition 14.11 We remark that this argument follows 
very closely with Proposition 2.9 in [6]. Nevertheless, we provide it here for completion. 

Proof, (of Proposition HH]) 1. Set 

£{r) := A^''^{r) + D'^ir), 0<r<l 

and observe 

rE{T) = r{{ -lin'] + r( -f4\P?'^^ ^ ^ -^^^ ^''^ ''^ 



<c{C{rf''' + D{rf/'^] 
<c{C{r) + D\r)Y'^ 

<c{C{r) + 8{r)Y'\ (4.6) 
Recall that our aim is to show that there is e > such that 

lim sup i?(r) < e ^ liminf ri?(r) < e^. 
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Therefore, we assume that hmsup^_j.o+ B{r) < e and choose e below to estabhsh the above 
imphcation. 

2. Let 6,p ^ (0, 1/2) be fixed and observe that easy corollaries of Lemmas 14.21 14.31 and 
14.51 are 

[lop) < c {c{ep) + D\ep) + A^/\ep)B^/^(ep)} , 
^D{ep) < c {eD{p) + e-^A^/\p)B^/\p)} . 

Using the inequalities (that are trivial to verify) 

A{9p)<^Aip) and B{9p)<^Bip), 

we also have 

A'^' (^^p) < c l^e'D^p) + e'A'/^p) + l-,Anp)Bnp) + + ^) Anp)B'/'{p) 



and 



( bp 



<c\e'D\p) + l-A''\p)B''\p) 



9^' 



With the above estimates on A^/"^ and -D^, 

s{^-ep)^A^i^{^-ep)^D^{^-ep 

< c [e^D^p) + e^A^I\p) + l-A^'\p)B^/\p) + 1a3/2(p)53/2(p) 



<c\e'£{p) + l,A'"{p)Bnp) + 

<<^[d'm+j,Bnp)] 

as A{p) < A{l/2) < c. 

3. Now choose 6' G (0, 1/2) so small that for each of the (universal) constants above 

c^^ < 1/2. 

We may also select 5^ G (0, 1/2) so small that 

B{p) < 2e for < p < 4 



(4.7) 



(4.8) 



25 



This can be done by hypothesis. With 6 G (0, 1/2) chosen to satisfy (14. Sp . (14. 7p gives 



S (rp) < ]^£{p) + Cr€"\ < p < 4 



where we have set 



r:= \e 



:=c23/V^9 ■ 

We interpret (14. 9 p to be a decay estimate for £^ as a simple induction argument provides the 
sequence of inequahties 

k 

^(^V) < ^Ap) + -re"' Y.h< I^Ap) + 2c.e^/^ (4.9) 

for any fixed p G (0, 5^) and G N. 
4. Employing 112] and (H^ . 

C(r'^+V) < c {r3A='/2(rV) + r-M3/^(rV)fi'/'(rV)} 
< c {r3£'(rV) + r-3£V2^rV)(2e)=^/^} 

'l£:(p) + 2c.e^/2) 



< c<^ r 



_-3 



2fc 



^(p) + 2c^€"' 



1/2 



,3/4 



In view of the above estimate and inequality (14. 6p . 

lim sup rV^(T^p) < dim sup [C(rV) +£{t^p)\^^^ 

< c |r^2c^ + r"^ (2c^)^/^ + 20^^^^ e 
We conclude by choosing e > so small that 

lim sup r''p£'(r''p) < ^-€5. 

fc— >-co ^ 



1/3 



J/2 



□ 



Remark 4.6. A close inspection of the work in this paper shows that the methods employed 
establish an analogous version of Theorem 11.11 for weak solutions of the system 



' dtu + (n ■ V)n = Au-Vp + V ■ FF* 
dtF + {u ■ V)F = AF + WuF + (VQ)* 
V ■ n = 
V ■ F* = 



(4.10) 
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where Q = Q{x,t) G 
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